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ABSTRACT 


This  report  presents  the  results  of  calculating  the  deflection 
and  stresses  in  a  thermally  loaded  beam.  Three  different  methods 
were  used  for  the  calculation:  classical  beam  theory,  ANSYS 
(Version  4.4)  and  VAASEL  (  Version  1.1) .  Comparisons  between  each 
of  the  methods  are  presented. 
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1 .  INTRODUCTION 


VAASEL  (Vulnerability  Analysis  of  Aerospace  Structure  Exposed 
to  Lasers)  is  a  finite  element  analysis  computer  code  developed  by 
the  Northrop  Corporation  for  the  Air  Force.  VAASEL  has  a  modular 
architecture  that  makes  it  a  valuable  tool  for  multi-disciplinary 
analysis.  VAASEL  is  currently  being  used  by  SAIC  (Science 
Applications  international  Corporation)  for  the  SAVE  (Structural 
Assessment  and  Vulnerability  Evaluation)  program. 

The  subject  comparison  came  about  when  SAIC  decided  to  run 
some  test  cases  through  VAASEL  and  repeat  these  cases  through  their 
current  techniques  for  laser  vulnerability  analysis.  SAIC  uses  the 
finite  element  program  ANSYS  to  do  the  thermal  and  structural 
analysis.  One  of  the  problems  chosen  was  a  trapezoidal  wingbox 
analysis  that  was  used  to  benchmark  the  VAASEL  code.  This  model 
contained  CSHEAR  elements  to  model  the  spar  webs.  Since  the 
current  version  of  ANSYS  does  not  support  shear  panel  elements , 
SAIC  replaced  the  shear  elements  with  membrane  elements  and  ram 
VAASEL  over.  This  resulted  in  a  maximum  tip  deflection  change  from 
+11.0  inches  to  +3.3  inches  for  a  combined  thermal  and  flight  load 
condition  (pure  thermal  loading  resulted  in  a  maximum  deflection  of 
-0.05  inches.)  This  change  in  deflection  was  expected  since  the 
CQDMEM1  elements  stiffened  the  model.  ANSYS  predicted  -6.6  inches 
deflection.  Although  it  was  not  independently  verify  that  SAIC  was 
running  the  same  problem,  it  was  recommended  that  additional 
modeling  be  done  to  check  the  validity  of  the  thermal  stress 
calculations  in  both  VAASEL  and  ANSYS  [1]  .  A  comparison  between 
the  two  codes  for  a  problem  with  a  classical  solution  was 
necessary . 


1 


2 .  ANALYSIS 


2 . 1  Problem  statement 

The  selected  problem  was  a  pure  bending  problem  with  a  beam 
theory  solution.  Figure  1  shows  a  schematic  of  the  idealized 
composite  beam.  The  beam  is  composed  of  one  material  but  each 
section  is  at  a  different  temperature.  The  reference  temperature 
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Figure  1:  Idealized  Beam  for  Thermal 
Load  Check 


was  taken  to  be  0°F.  Under  this  loading,  the  top  and  bottom  truss 
beams  will  try  to  expand  while  the  web  of  the  beam  will  restrain 
them  causing  the  beam  to  bend. 

2.2  Classical  Solution 

A  strength-of-materials  analogy  is  used  to  determine  beam 
deflection  [2].  This  formulation  ignores  shear  stress  in  the  beam. 
The  closed  form  solution  is  shown  in  Appendix  A.  The  results  are 
listed  in  Tables  I-II. 


2.3  Finite  Element  Solutions 
2.3.1  VAASEL  Model 

To  check  VAASEL,  two  models  were  run:  one  containing  CQUAD4 
elements,  and  the  other  containing  QDMEM1  elements.  The  QUAD4 
element  is  an  isoparametric  quadrilateral  plate  element  with  both 
membrane  and  bending  behavior.  The  QDMEM1  is  an  isoparametric 


quadrilateral  membrane  element . 
Both  elements  were  tested  since 
SAIC  used  the  QDMEM1  elements  for 
the  trapezoidal  wing  box  problem, 
but,  the  element  used  for  the 
ANSYS  analysis  is  more  closely 
related  to  the  QUAD 4  element. 
Both  VAASEL  models  consisted  of  1 
row  of  CROD  elements  along  the  top 
and  bottom  edges  to  model  the 
truss  beam  (Figure  1)  .  The 
membrane  web  was  modelled  with  two 

rows  of  quadrilateral  elements . 

The  current  version  of  VAASEL 
doesn't  support  element 
temperature  assignment  (this 
problem  was  discovered  while  doing 
the  analysis  and  will  be  corrected 
in  the  next  version  of  VAASEL)  . 

To  get  around  this  problem  all 
temperatures  were  assigned  to  the 
grid  points  and  the  membrane  web 
was  given  a  0.0  coefficient  of 
thermal  expansion.  This  forced 
the  solution  to  look  like  the 
classical  problem  (no  thermal 
expansion  in  the  membrane  web) . 

2.3.2  ANSY8  Model 

The  ANSYS  model  was  run  by 
Mr  Chuck  Heighland  of  SAIC  [1]  . 

His  model  used  2  rows  of  STIF63 
(quadrilateral  shell)  elements  to 
model  the  membrane  web  and  a  row 

of  STIF8  (spar)  elements  along  the  TABLE  II:  Displacement 
top  and  bottom  to  model  the  truss  Results 

beam.  ANSYS  supports  element 

temperature  assignment  and  this  was  used  to  analyze  the  beam. 


Table  I :  Stress 
Results 
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3 .  RESULTS 


The  VAASEL  and  ANSYS  results  compare  very  well,  all  results 
are  within  2%  (Table  III).  The  VAASEL' s  models  were  slightly 
stiffer  particularly  with  the  CQDMEM1  card.  The  additional 
stiffness  using  the  CQDMEM1  card  was  predicted.  The  character  of 
the  finite  element  solutions  were  as  expected.  The  axial  stress  in 
each  component  was  nearly  const  amt,  except  for  end  effects,  in 
addition,  the  axial  displacements  were  linear  and  the  vertical 
displacements  were  quadratic  (Figure  2  &  3)  .  The  comparison  to 
beam  theory  appears  to  be  good,  since  beam  theory  disregards  shear 
stress,  while  both  ANSYS  and  VAASEL  include  all  stresses.  Appendix 
B  contains  all  the  displacement  results. 


Description 

Beam 

Theory 

ANSYS 

VAASEL 

(CQUAD) 

VAASEL 

(CQDMEM1) 

Ux  (inch) 

.040 

.033 

.033 

.033 

Uy  (inch) 

.834 

.  654 

.645 

.552 

Oi  (pai) 

-6167 

-5844 

-5736 

-5288 

Ui(pai) 

-9832 

-11298 

-11284 

-11136 

<j3(psi) 

4000 

3429 

3422 

3423 

TABLE  IZZ:  Summary  of 
Results 
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APPENDIX  X 


This  section  describes  the  classical  solution  to  the  thermally 
toaded  beam  problem.  The  primary  equations  used  are  listed. 


JaxzdA=bjoxzdz 
P=  joxdA=bjaxdz 


Ox-Ei(u!0-Z'Hn  -Ota  T±) 


The  analysis  and  results  are  shown  below. 


—  •  JaxzdA=bjoxzdz-bjEi  ( u7  -zw"  -a1T1)  zdz 


Figure  X-l :  Problem 
Setup 
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Let  AT,  *  0  and  E,  ■  E,  ■  B, 
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A-l 


Integrate : 


8  n 


u‘  ~A 
u  =Ax+C 


Apply  end  constraints: 


u=0  @x=0  =>  C=  0 


u=Ax 


(ATx+AT3) 


1-. 


i) 


Calculate  the  deflections  and  stress  using  the  following  numbers 

-  10  Mai  ho  -  2.0  Tx  -  100°F 

-  6  E-6  in/in- °F  h3  -  3.0  T3  «  300°F 


X  (in) 

u*(in) 

Uy(in) 

10 

0.004 

0.008 

20 

0.008 

0.033 

30 

0.012 

0.075 

40 

0.016 

0.133 

50 

0.020 

0.208 

60 

0.024 

0.300 

70 

0.028 

0.408 

80 

0.032 

0.533 

90 

0.036 

0.675 

100 

0.040 

0.834 

Table  A-l:  Displacement  Results 


.-3 


<*i  (pai) 

-6167.5 

<y2  (psi) 

-9832.5 

a3  (psi) 

4000 

Table  A-2:  Stress  Results 
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APPENDIX  B 


DISPLACEMENT  RESULTS 


x  (in) 

u,  (in) 

u,  (in) 

x  (in) 

us  (in) 

WBM 

0 

0.000 

0.000 

0 

0.000 

0.000 

10 

0.003 

0.006 

10 

0.003 

0.005 

20 

0.007 

0.026 

20 

0.007 

0.022 

0.010 

0.058 

30 

0.010 

0.050 

40 

0.014 

0.103 

40 

0.014 

0.088 

50 

0.017 

0.161 

50 

0.017 

0.138 

60 

0.020 

0.232 

60 

0.020 

0.198 

70 

0.024 

0.316 

70 

0.024 

0.27 

80 

0.027 

0.413 

80 

0.027 

0.353 

90 

0.031 

0.522 

90 

0.031 

0.447 

100 

0.033 

0.645 

100 

0.033 

0.552 

Table  B-l :  QUAD 4 
Displacement  Results 


Table  B-2:  CQDMEM 
Displacement  ^ Results 
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